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Exploiting the recently derived quantum Smoluchowski equation the classical 
Ivanchenko Zil'berman theory for overdamped diffusive phase motion of low ca- 
pacitance Josephson junctions is extended to the low temperature quantum domain 
where charging effects appear. This formulation allows to derive explicit results for 
the current-voltage characteristics over a broad range of parameters that reduce to 
known findings in certain limits. In particular, the transparent analytical approach 
comprises Coulomb blockade physics, coherent Cooper pair transfer, and the pre- 
cursors of macroscopic quantum tunneling and needs to be supplemented by more 
sophisticated methods only at very low temperatures. 



I. INTRODUCTION 

In the last decades the physics of Josephson junctions ( JJ) has revea! 
wealth of phenomena studied theoretically and experimentally as well 
lying system, two superconducting domains separated by a tunnel barrier, can be found in 

n 

a variety of realizations, recently e.g. in superconducting atomic contacts •4] and solid-state 



ed an extraordinary 
l|,3,l3. The under- 



quantum bits 5]. Basically two parameters determine the dynamics, namely, the Josephson 
coupling energy Ej and the charging energy = 2e^/C of a junction with capacitance 
C . It turns out that the competition between these two scales is crucially influenced by 
the electromagnetic environment surrounding the junction, i.e. its impedance which in the 
simplest case is given by an ohmic resistor with resistance R. 

For the classical JJ dynamics charging effects are negligible and the charge transfer cor- 
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responds to a phase- coherent Cooper pair current 

/s = /esin(0) , 0(t) = ^T/(t) (1) 

with Jc = {2e/h)Ej, the phase difference across the junction 0, and where the phase velocity 
(■ = d/dt) is related to the voltage drop V{t) across the junction. In the opposite limit of 
dominating Ec, however, charging effects prevail and Coulomb blockade is known to rule the 
transfer of incoherent tunneling of Cooper pairs j^. 

A particular appealing feature of a JJ is the fact that its classical dynamics can be visu- 
alized as the diffusive motion of a fictitious classical particle (RSJ model) Accordingly, 
for a current biased JJ the case of very small capacitance corresponds to very strong friction 
(Smoluchowski limit) associated with the phase dynamics 

m-f^ + dU{(f))/d(f) = ^{t) . (2) 

Here the translation rules are: Mass m = (/i/2e)^C, damping strength 7 = 1/RC, potential 
U{(j)) = —Ej cos{(f)) — E]^(f) with the energy E]^ = {h/2e)I related to the bias current, and 
current noise (^(t)) = 0, {^(t)C,(t')) = (27//?) 5{t — t') {[3 = l/k^T). The corresponding time 
evolution equation for thephase distribution (Smoluchowski equation) has been studied in 
detail already in the 60s jy, 3]- Even in the overdamped limit the charging energy E^ can 



ne 



be relatively large and the crucial question arises in which way charging effects alter t 
diffusive phase dynamics of JJs. While this has been analyzed at very low temperatures 
and for small and large [l^ Josephson couplings, much less is known for intermediate 
temperatures and coupling energies. 

The goal of this paper is to fill this gap for overdamped junctions by providing a gen- 

n 

eralization of the classical Ivanchenko-Zil'berman Theory (IZT) 6]. As we will show, the 
classical and quantum regimes are characterized by /SEc/np <^ 1 and /SEc/np ^ 1, respec- 
tively, with p = R/Rq {Rq = h/Ae^) meaning that in a mechanical analog where 

f5Ej'Kp = n(3-f (3) 

the appearance of charging effects is associated with the changeover from classical (7/1/^ -C 1) 
to quantum (7^/? ^ 1) overdamped phase diffusion. While in general the dissipative dy- 
namics of a quantum particle is a formidable task ^|, it was shown recently [12] that in 
the overdamped limit considered here the reduced quantum dynamics can be described by 
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a general extension of the classical Smoluchowski equation, the so-called Quantum Smolu- 
chowski Equation (QSE). This allows to reveal the relation between classical and Coulomb 
blocked transport in overdamped JJs in terms of a simple analytical approach. 



II. QUANTUM PHASE DIFFUSION IN JJS AT STRONG FRICTION 



Classical overdamped motion is based upon a separation of time scales between relaxation 
of position and relaxation of momentum Q]. Momentum is slaved to position and on a 
coarse grained time scale can assumed to be in thermal equilibrium with respect to the 
instantaneous position of the particle. 

Quantum dissipative dynamics is much more complicated and usually has to be formu- 

T I 

lated in terms of the path integral formalism [11]. This approach has been exploited very 
successfully in the past to describe low temperature propert.es of JJs in parfcular 030 

and of dissipative particles in periodic potentials in general [HI, Explicit calculations 
require involved techniques though, and analytic expressions are available only in certain 
ranges of parameter space. The problem is that quantum fluctuations acting on time scales 
of order hjS lead, particularly at low temperatures, to long range interactions in time jll |. 
Recently it was shown, however, that in the overdamped domain the exact quantum time 
evolution can be mapped onto a quantum generalization of the classical Smoluchowski equa- 
tion for the position probability distribution. The resulting QSE reads 



dP{4>,t) 



1 d 



(4) 



dt 7m d(j) 

with ' = d/d(j) and the translation rules as specified above. Quantum fluctuations appear as 
a Josephson potential with an effective Josephson coupling 



[/eff(0) = -^j* COS(0) - Eb0 , = E, 



1-A 
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and a phase dependent diffusion coefficient capturing quantum noise jl^ 

D{(j)) = [1 - e cos{(f))]-\ e = A(3Ej. 



(5) 



(6) 



The changeover from classical to quantum dynamics is governed by a function A that contains 
'jhp as the essential quantity. With (jH} it takes the explicit form 



A = 2p 



27c^p ^ ( pE, 
c+— f + ^ ' 



I3E, 



27r2p, 



(7) 



-7r/2 TT Itt 



FIG. 1: Effective potential Ues{(l>) for Ey, = 0.4 (solid, in units of ii^j) and diffusion coefficient D{(f)) 
(dashed, for (3Ej = 1) together with the bare Josephson potential U{(j)) (dotted, in units of Ej). 
Other parameters are PEc = 20, p = 0.04. 

where c = 0.5772 ... is Euler's constant and \l/(-) denotes the logarithmic derivative of the 
gamma function. In the high temperature limit l3Ec/ p <^ 1 quantum fluctuations become 
small and independent of p, i.e., A ^ f3Ec/'n''^, and one regains classical Smoluchowski 
dynamics. In the low temperature range (3Ec/ p ^ 1 one has A ^ 2plog(/9i?c/27r^p) leading 
to a substantial influence of quantum fluctuations. Qualitatively, the effects of -Ej and 9 on 
the dynamics of the fictitious particle are easily understood (cf. fig. QJ: E'j < i?j lowers the 
potential barrier in Ucs-, ^-dependent terms in the diffusion coefficient [see (jSEI)] enhance the 
noise strength in the wells {D > 1) and lower it around the barrier tops {D < 1). It turns 
out that this causes a moving particle to experience a larger mobility where in limiting cases 
barrier (/J-Ej ^ 1) or diffusion (/5i?j ^ 1) related quantum fluctuations dominate. 

Let us now first address the range of applicability of ^ for the physics of JJs. The 
QSE is valid if (i) a separation of time scales hjS,!/'-/ <^ l/^j with the plasma frequency 
tuj = {2e/h)RIc is guaranteed. Classically, for the McCumber parameter this condition 
reduces to jSj = j /ujj ^ 1. Further, (ii) in order for the momentum m</) to relax within the 
RC = 1/7-time to a Boltzmann distribution around (0), the external voltage V is restricted 
by eV <^ h'-f. We note in passing that (ii) also ensures that the actual non-ohmic impedance 
seen by the junction can effectively be treated as ohmic. By combining (i) and (ii) and 
re-expressing them in junction parameters we thus expect the QSE to describe quantum 
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phase diffusion in J J if (see fig. |2I) 



» 1, npPE,, ^ . (8) 



Ej27r2p2 ' RIc 
Since typically p <^ 1 the above condition allows for a broad range of values for Ec/Ej, (3Ej, 
and also large voltages V/RIc- It is important to note that for realistic J J with finite p the 
classical description 0, ^ is always restricted to small and moderate values of PEj, while 
with increasing one approaches the quantum Smoluchowski range with a finite A [see 
(jZj) and fig. El- 

To complete this discussion we also specify the corresponding quantum Langevin equation 
(in the Ito sense) which generalizes the classical one Q to 

m^^ + UM)=m v^W)- (9) 

Now, the Eqs. (jU and Q are the starting point to study in a very elegant way all aspects 
of the phase dynamics of JJs within the range (jH)). While in the quantum domain very low 
temperatures (T 0) have been explored in detail ^^J, finite temperature results exists only 
for small and large [l^ Josephson couplings, respectively. In the sequel we will show that 
the QSE bridges between these limiting cases. 



III. STEADY STATE CURRENT AND CURRENT VOLTAGE 

CHARACTERISTICS 



For finite bias current If, > the phase dynamics Q approaches a stationary non- 
equilibrium state at longer times. This is associated with a stationary current Jst = 
hmt_oo(0(t))/27r given by 



J. 



St 



1 _ q-P^-kE^, 

7m/3 



2tt 



Di<l>) 



with 



^{(P) = -Ej^ cos(0) - E^<p + OE^ sin(0) - sm\ct>) . 



(10) 



(11) 



According to {V) = (/i/2e) Jst one then finds the current voltage-characteristics of a current 
biased junction to read 



(V) 



pn 1 



(12) 
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FIG. 2: Range of the QSE for a JJ with p <C 1, V/RIc < 1- The classical range (shaded) and 
the domains of Coulomb blockade (CB) and macroscopic quantum tunneling (MQT) are indicated. 
The QSE is applicable above the thick line, see ((SI), and the arrows illustrate various changeovers 
discussed in the text. The slope of the CB boundary decreases with increasing V/RIc so that for 
V/RIc > 1 CB dominates in most parts of parameter space. 

where the nominator Tq^ results from normalizing the steady state phase distribution to 1 
and can be written as 



qm 



1 /-s^ 
27r io 



2tt 



g-2/3£J cos(</,') sin(0/2) _ q _ ^/2)] ^WeH^P,^') (^g) 



with 



/,/g-/3i?b<^'j^[2/?£;jsin(072)] 



(15) 



e(0, 0') = sin(0') sin(0/2) [E^ + E', cos(0') cos(0/2)] . (14) 

The expression (jl2p together with (jl3p is the central result and will be discussed in the 
following. 

For PEc/p <ti 1 the function Tqm reduces to the classical result 

/■2tt 

Tci{Ej) = / , 
Jo 

with /q the modified Bessel function [l| and IZT j6|, l2| is recovered. On the other hand, 
in the low temperature domain jSEc/ p ^ 1 two distinct ranges must be distinguished (see 
fig. El). For smaller couplings, jSEj < 1, where Coulomb blockade (CB) transfers charges 
incoherently, we have ^ 1 so that diffusion related quantum fluctuations are negligible. 
Then, fllHj] takes the classical form with quantum fluctuations incorporated in terms of a 
renormalized Josephson energy, i.e. Tq^ ~ Tci(-Ej). This important extension of IZT has 



first been derived in 0] based on a direct evaluation of tlie real-time path integral expression. 
However, since diffusion related quantum fluctuations are disregarded, this simple result is 
(for finite p <^ 1) indeed restricted to In contrast, (fT!^ also holds for much larger 

E] [fig. 121 arrow (a)] associated with coherent Cooper pair tunneling. An explicit formula 
can be obtained from ()12|1 for sufficiently large jSEj > 1 and a = Ei,/ Ej = I/Ic < 1. Then, 
occasional phase slips lead to the voltage 



(y)_ ^ yi - -2/3Ej(l-a2)3/V(3a2) 2eVT^ (IQ) 

RIc 2tt ^ ' 

which via 9 is affected by diffusion related quantum fluctuations. Accordingly, quantum 
effects cannot always be captured by a renormalization of the Josephson coupling, rather such 
a simplification occurs for CB dominated transport only. As can be observed in fig. |21 the 
result (fT^ tends for ^ to classical thermal activation over the barriers of the washboard 
)otential f/(0) where quantum corrections are of the known damping independent form 
IgI ]. At lower temperatures, i.e. for finite 6', they show a complicated dependence on p 
and capture the precursors of macroscopic quantum tunneling (MQT) found at very low 
temperatures Thus, the central result ()12|1 fills the gap between established results in 
different transport domains: On the one hand for fixed (3Ec/ p > 1 it leads with increasing 
from Coulomb blockade to coherent Cooper pair tunneling [fig. El arrow (a)]. On the 
other hand for fixed > 1 it connects with varying (3Ec/ p classical thermal activation 
with MQT [fig. arrow (b)]. Apart from limiting cases, the result (fT^ must be evaluated 
numerically as illustrated in fig. IHl (left). 

What is left, is the domain of small Josephson coupling < 1 and increasing jSEc/p 

for a voltage biased junction. 
17[ . Here, we demonstrate that 
the QSE recovers this scenario, too, and hence gives a complete description throughout the 
QSE range sketched in fig. EJ Results for a voltage biased JJ are gained from the current 
biased case by means of / = V/R and {V) = V — RIs- The interesting observable then is 
the supercurrent Is- From ()12|1 we obtain to order (/?i?j)^ 

Je 2 Rhl + [PE,V/{Rh)f ^ ^ 

with the characteristic renormalized E;] discussed above. This result coincides with that 
of Q and connects quite different transport regimes [see figs. I2|3I (right)]. Namely, in the 



[fig. El arrow (c)l that has been discussed in the literature la 

n 

This changeover leads from classical coherent IZT 6] to CB 
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FIG. 3: Left: Current-voltage characteristics for PEj = 2, p = 0.04. The quantum results [PE(. = 1 
(dashed), PE^ = 20 (solid)] are shown with the classical one {(3Ec = 0, dotted). Right: Supercur- 
rent vs. voltage for = 0.25, p = 0.04. The classical result (dotted, j^) and the CB expression 
[dashed, (0] are depicted together with the full QSE resuh (sohd) for l3Ec = 0.15 (top), (3Ec = 20 
(middle), PEc = 1000 (bottom). 



classical limit (3Ec/ p <^ 1 it gives the small-/? i?j-limit of the IZT while in the opposite 

H.ntum domain l3Ec/ p ^ 1 it agrees with the small p, A expansion of the CB expression 

h f3Ej e"^ \r{p-tPeV/7r)\' . 

— = z — ^ smh(BeV) . (18) 

/c 47r r(2p) ^ ^ ^ 

Thus, we conclude that the quantum Smoluchowski result (fTTjl interpolates in the domain 
of small Josephson couplings between phase- coherent and incoherent charge transport (with 
E] Ej) where the strength of charging effects controls the changeover from one regime 
to the other (see fig. Note that p7|) yields the supercurrent reasonably well also for 
/3Ec > 1 (A of order 1) [fig. H (right), bottom]. 



IV. CONCLUSIONS 



Based on the analogy between JJ physics and the diffusive motion of a fictitious particle in 
a washboard-type of potential we generalized for strong friction (low capacitance) the IZT to 
the low temperature quantum domain. Quantum fluctuations affect the potential profile as 
well as the diffusion constant and physically are associated with charging phenomena. This 
formulation reduces to known results derived in certain limits and thus gives a complete 
description over broad domains of parameter space. Our main expression is specified in fjl2|) 
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together with ()13p. It contains incoherent charge transfer (CB) and coherent Cooper pair 
tunnehng as well as the precursors of MQT. In particular, the incorporation of quantum 
fluctuations in terms of a renormalized Josephson junction turned out to be characteristic 
for CB only. Eventually, our findings confirm the recently developed quantum Smoluchowski 
theory to be a powerful tool for studying low temperature quantum diffusion in mesoscopic 
physics. 
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